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Aims. We present an analysis of OGLE 2004-BLG-254, a high-magnification (A a 60) and relatively short duration (/e — 13.2 days) 
microlensing event in which the source star, a Bulge K-giant, has been spatially resolved by a point-like lens. We seek to determine the lens 
and source distance, and provide a measurement of the linear limb-darkening coefficients of the source star in the / and R bands. We discuss 
^ C~| the derived values of the latter and compare them to the classical theoretical laws, and furthermore examine the cases of already published 
O i 1 microlensed GK-giants limb-darkening measurements. 

J- Methods. We have obtained dense photometric coverage of the event light curve with OGLE and PLANET telescopes, as well as a high 
t| . signal-to-noise ratio spectrum taken while the source was still magnified by A ~ 20, using the UVES/VLT spectrograph. We have performed a 
' modelling of the light curve, including finite source and parallax effects, and have combined spectroscopic and photometric analysis to infer 
the source distance. A Galactic model for the mass and velocity distribution of the stars has been used to estimate the lens distance. 
JL* Results. From the spectrum analysis and calibrated color-magnitude of the event target, we found that the source was a K3 III Bulge giant, 
situated at the far end of the Bulge. From modelling the light curve, we have derived an angular size of the Einstein ring d E = 114 ^as, and a 
relative lens-source proper motion // = #eAe — 3.1 mas/yr. We could also measure the angular size of the source, 0, m 4.5yuas, whereas given 
5^ ■ the short duration of the event, no significant constraint could be obtained from parallax effects. A Galactic model based on the modelling of 
Ci ' the light curve then provides us with an estimate of the lens distance, mass and velocity as Z\ - 9.6 kpc, M at 0.11 M Q and v - 145 kms -1 
(at the lens distance) respectively. Our dense coverage of this event allows us to measure limb darkening of the source star in the / and R 
bands. We also compare previous measurements of linear limb-darkening coefficients involving GK-giant stars with predictions from ATLAS 
atmosphere models. We discuss the case of K-giants and find a disagreement between limb-darkening measurements and model predictions, 
which may be caused by the inadequacy of the linear limb-darkening law. 
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1. Introduction 

The microlensing technique is one of a few methods (together 
with interferometry, eclipsing binaries and transiting extra- 
solar planets) that can be used to measure brightness profiles 
and limb-darkening coefficients of stars at distances exceed- 
ing a few kpc. The apparent stellar disk is a projection of the 
near stellar hemisphere and therefore shows an axisymmetric 
variation of the observed brightness as the projection maps 
the variation of the emergent radiation with depth. The con- 
tinuum is on average formed at larger depth at the disk cen- 



ter and a t small er dept h at the limb (Eddington-Barbier effect, 
see e . g.: IWitt & Maolll994l " 
ll997tlHevrovskv et al.l200(X 



Loeb & Sasseloylll995t Isasselovl 
Hevrovskvl2003l) . The microlens- 



* Based on observations made at ESO, 073.D-0575A 



ing method requires the source to transit a region with a large 
magnification gradient over its face, as present in the vicinity 
of caustics. While a single lens creates a point-like caustic at its 
angular position, binary lenses create extended caustic patterns 
formed of lines (fold caustics) which merge at cusps. 

PLANET observations of the m i crolensi ng event 
MACHO 1997-BLG-28 jAlbrow et al] Il999bl) contain- 
ing a cusp passage, constituted the first limb-darkening 
measurement of a Bulge giant. The majority of subsequent 
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limb-darkening measurements such as MAC HO 1998-SMC-l 
jAlbrow et all Il999al lAfonso et alJ Eoool MACHO 1997- 



(lAlbrow et al.l 


2001 


), EROS-2000-BLG-5 
, OGLE 2002-BLG-069 


2003 


:lAnet al.ll2002l 


2004 


: IKubas et 


al l 120051 iThurl et al.l 12006 ) 



fold-caustic passages. The limb-darkening measurement in 
the solar-like star MOA 2002-BLG-33 JAbe et alJ l2003h 
constituted a very special case with the source enclosing 
several cusps of the caustic at the same time. In the single lens 
case, the extended size of the source star has a significant effect 
on the light curve if the angular source size is of the order 
or larger than the angular separation between source center 
and lens point-like caustic. However, only a small fraction 
of microlensing events provide angular separations that are 
small enough , and so far very few cases have been obs erved 
llAlcocketal.lll997t lYoo etaljEooH I.Tianget aljEool) . For 
events showing evidence of this effect, the limb darkening of 
the source star can be determined, and constraints on the phys- 
ical properties of the lens can be derived using estim ates of 
its dis t ance, spectr a l type and event model parameters JGouk 
1994 [ IWittl Il995t lLoeb & Sasselovl Il995t | Hevrovskv etal 



2000: lHevrovsk 



v 2003). 



This work treats the case of OGLE 2004-BLG-254, a high 
magnification event from a point-like lens showing extended 
source effects, for which a dense photometric follow-up was 
performed at four PLANET observing sites, making it one 
of the best observed events of this kind to date. Using high- 
resolution spectra collected with UVES (VLT), we derive the 
characteristics of the source star, a K3 giant. We search for a 
suitable single-lens limb-darkened source microlensing model, 
which we use together with the properties derived from the 
high resolution spectra to discuss the source center-to-limb 
variations and constraints on the lens properties. We compare 
our limb-darkening measurements to previous observations of 
such stars, especially to EROS-2000-BLG-5 for which a large 
discrepancy was found between observations and atmosphere 
models. We also note that the derived characteristics provide 
the basic information necessary for a forthcoming study of ele- 
ment abundance in the Bulge giant source star of OGLE 2004- 
BLG-254. 



2. Photometric and spectroscopic observations 

2.1. Photometric monitoring 

The OGLE-III Early Warning System (EWS) JUdalskill20oi 
discovered and alerted the Bulge event OGLE 2004-BLG- 
254 (a = 17h56m36.20s, 5 = -32°33'01'/8 (J2000.0) and 
I = 358.09°, b = -3.87°) on May 17, 2004, from observa- 
tions carried out with the 1.3 m Warsaw Telescope at the Las 
Campanas Observatory (Chile). 

The PLANET collaboration started its photometric ob- 
servations on June 8 which form the basis for our anal- 
ysis and consist of data from 4 different telescopes being 
part of the PLANET network: the Danish 1.54m at ESO La 
Silla (Chile), the Canopus lm near Hobart (Tasmania), the 
Elizabeth lm at the South African Astronomical Observatory 



(SAAO) at Sutherland and the Rockefeller 1.5m of the Boyden 
Observatory at Bloemfontein (both in South Africa). Every 30 
minutes the data from the different sites are uploaded to a cen- 
tral computer in Paris, where data are combined and fitted au- 
tomatically, and light curves are made publicly available. 

The event was also monitored by the //FUN collaboration 
from Chile (1.3m telescope at the Cerro Tololo Inter- American 
Observatory) and Israel (Wise lm telescope at Mitzpe Ramon). 
Although we did not make use of them for our modelling, we 
checked they were consistent with our data sets. 

The data we collected on OGLE 2004-BLG-254 showed a 
rise in apparent brightness by 2.85 mag above baseline on June 

9, 8:10 UT. These pre-peak data and adequate modelling pre- 
dicted a peak to occur on June 10, 6:35^3g™ UT, at a rather 
uncertain, but in any case large, magnification of A a = 80^]. 
Events of this type harbour an exceptional potential for the dis- 
covery or exclusion of extra-solar planets as well as for the 
study of stellar atmospheres and might provide an opportunity 
for measuring the mass of the lens star. 

On June 10 at 12:45 UT, a public alert was issued by 
PLANET reporting that data collected on OGLE 2004-BLG- 
254 at the SAAO 1.0m between June 9, 18:50 UT and June 

10, 4:40 UT, at the Danish 1.54m on June 10 between 2:20 UT 
and 10:05 UT, as well as OGLE data obtained on June 10 be- 
tween 3:50 UT and 9:55 UT, revealed the extended size of the 
source star, and the crossing time of the disk was estimated to 
be around 16 hours. The peak was passed around June 10, 7:40 
UT at 4.35 mag above baseline, at a magnification of A ~ 55. 



2.2. Blending of the source star 

There are two nearby stars on the eastern side of OGLE 2004- 
BLG-254. Their calibrated OGLE magnitudes and colours are: 
/ = 16.31, (V - I) = 2.30 for the star located 0.23" north and 
1.41" east, and / = 16.69, (V - I) - 1.55 for the star located 
0.52" north and 3.00" east. 

Thanks to its relative brightness, the source star of 
OGLE 2004-BLG-254 c an be found in rec ent infrared 
surveys, such as DE NIS fepchtein et ail 1 1994 or 2MASS 
JSkrutskie et alJll997h . However, due to the large pixel size of 
these surveys (1 to 3") and the proximity of a companion of 
similar colour, the infrared measurements probably correspond 
to a blend of these stars: DENIS measured a "star" at 1.0" E 
and 1.5" S from OGLE 2004-BLG-254, with / = 15.81±0.09, 
J = 14.02 ± 0.12, K s = 13.12 + 0.16, while 2MASS measured 
the same "star" at 1.3" E and 1.0" S, with J = 14.135 ± 0.046, 
H = 13.326 + 0.048 and K s = 12.969 ± 0.045. 2MASS quality 
flags are optimal, and blend and confusion flags are not acti- 
vated. The DENIS / magnitude is therefore 0.78 mag brighter 
than the OGLE calibrated measurement (cf. Sec. 14. 1 1 , which 
supports our hypothesis that DENIS measured a blend of the 
two nearby stars of similar magnitudes. 



2.3. UVES spectroscopy 

On June 11, 2004, between 00:24 and 00:52 UT (HJD 
= 2,453,167.54159 at mid-exposure), we obtained a high- 
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resolution spectrum of OGLE 2004-BLG-254 using the U VES 
spectrograph mounted at the Nasmyth focus of Kueyen (VLT 
second unit). At this time, the source flux was magnified by the 
gravitational lens by a factor 20, making the VLT equivalent to 
a ~ 37 m diameter telescope. The spectrum was taken at one 
of the standard red settings centered on 5800 A, which cov- 
ers the spectral domain 4780 - 6835 A at a mean resolution of 
~ 40, 000, and typical S/N of 100. 

The data were reduced using version 2.1 of the UVES con- 
text of the MIDAS data reduction software. The raw science 
data were first bias-substracted and then wavelength-calibrated 
order by order using a Thorium-Argon lamp. The position of 
the science spectrum was determined in each order and an op- 
timum extraction was used to obtain ID science spectra where 
the sky lines had been removed. An Halogen lamp was used to 
obtain flat-field spectra, which were bias-substracted and com- 
bined, then wavelength-calibrated. The spectra were then flat 
fielded order by order, and merged to obtain the final data. The 
detail analysis of the final spectrum is detailed in Sec. 14.21 

3. Modelling of the light curve 

The photometric data collected on OGLE 2004-BLG-254 (see 
Fig.[D clearly show that the light curve is affected by extended 
source effects. 



3.1. Extended-source formalism 

With Ds denoting the source distance and x = D^/Ds the frac- 
tional distance of the lens, of mass M, the angular Einstein ra- 
dius reads 



4GM 1 -x 



(1) 



which is a characteristic scale of microlensing. For a point 
source situated at a projected angular separation u 0% from the 
lens, the magnification function is given by 

u 2 + 2 

Apspl(m) = ■ (2) 

u Vm 2 +4 

The magnification Aes of an extended limb-darkened source of 
angular radius p, 6>e is obtained by integrating Apspl over the 
source disk. With I{r) being the brightness profile of the source, 
normalized to unit flux (Sect. l3~2l . 



A E s(u\p*,I) 



where 



j j I(r) A PSPL (?7) r dtp dr, 
o o 



(3) 



"> 2 - 2 u rp„ cos tp + u 2 . 

For a uniformly bright source, IWitt & Mad 1 1994 have de- 
rived a semi-analytic expression of Aes, involving elliptic inte- 
grals. However, based on the fact that in a usual microlensing 
event toward the bulge of the Milky Way, extended source ef- 
fects are only prominent for small le ns-source angu lar separa- 
tions (« <K 1), where Apspl(m) — u -'. lGouldNl994l) found that 
the extended-source magnification factorizes as: 



(4) 
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Fig. 1. The upper panel shows the photometry of the mi- 
crolensing event OGLE 2004-BLG-254 near its peak, on June 
10, 2004, observed by four PLANET sites, Danish 1.54m, 
UTas lm, Perth 0.6m, Boyden 1.5m and SAAO 1.0m and 
OGLE. Since the shape of the magnification curve depends 
on the value of each linear limb-darkening coeffecient, differ- 
ent for each telescope, we choose the common best value of 
UTas and OGLE in Tab. ^ t° plot (solid line) the best-fitting 
point-lens, limb-darkened extended source model around the 
peak region. The residuals for each site are displayed by as- 
suming their own best linear limb-darkening coefficients. The 
lower panel shows the residuals of the complete set of data (the 
two vertical lines indicate the peak region displayed above). 
The UVES spectrum was taken at MHJD= 3 167.54159 at mid- 
exposure (vertical arrow). 



The second factor can be expressed by the s emi-analytical fo r- 
mula B Q (z) = iz£[arcsinmin(l,i),z] JYoo et all EooH) , 
where E is the incomplete elliptic integral of the second kind. 
By separating the u and z = u/p* parameters, this formula 
allows easy discretization and fast computation of extended- 
source effects. 

For non-uniform profiles (e.g. power-law limb-darkening 
models or tabulated profiles from stellar atmospheric mod- 
els), no such simple exp ressions are known. Different strate- 
gies have been proposed: IWitt & Maol d 19941) use the uniform 
source magnification and its derivative in a one-dimensional 
integral; iHevrovskvl d2003l) first calculates analytically the an- 
gular integral of the magnification, so that a single integral in- 
volving the (ra dially dependent) brightness of the source re- 
mains. Finally, lYoo et alJ d2004 give expressions of B r and 
Bi/2- functions, related to the linear and square root limb- 
darkening laws, respectively (with the same approximations as 
for Eq. to be numerically integrated. 
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When applying this formalism to OGLE 2004-BLG-254, 
we find that the maximal relative discrepancy between the ex- 
act magnification and its approximation related to Eq. (0} is 
less than 0.05 % for a linear limb-darkening law, which is well 
within the typical data error bars. We therefore use Eq. d4} to 
derive our limb-darkening coefficients. 

3.2. Event model parameters 

In the photometric analysis of the event, both PLANET and 
OGLE data are used. For each PLANET observation site, we 
have applied a cut on seeing which only removes very unreli- 
able points; we restricted the complete OGLE data to the ones 
collected after HJD' = 3050.0 (which is large enough to de- 
rive the baseline magnitude). OGLE provides 128 data points, 
SAAO 1 14, UTas 64, Boyden 75 and Danish 231, for a total of 
5 observing sites and 612 measurements. Data reduced with 
our DoPhot-based pipeline underestimate photometric errors 
(e.g. for bright magnitudes, error estimates may reach 10~ 3 , 
which is unrealistic). Comparison with the scatter of the non- 
variable stars suggests that the errors are underestimated by 
about 20%. Based on PLANET experience of PSF photometry 
used here, we rescale the errors as <x 2 esc = (1.2 cr data ) 2 + (0.01 ) 2 . 
As the source star is a K3 giant, we first check if the fluctu- 
ation in the OGLE baseline magnitude can be due to intrin- 
sic periodic activity (which could easily be included into the 
models). A power spectrum of 115 baseline data points com- 
ing from OGLE does not show evidence for modulation greater 
than cr mod = 0.015. We therefore attribute the residuals in the 
baseline to non-periodic variability or observational noise. 

As expected, a single-lens point-source model results in a 
high^ 2 /d.o.f. (~ 31), which obviously excludes this possibil- 
ity. On the other hand, a uniform extended source model pro- 
vides us with a first working set of parameters that fit the data 
well. However, at this stage, the residuals of the fit show some 
symmetric trends around the peak of the light curve that indi- 
cate limb-darkening effects. We thus add linear limb darkening 
to the source model. This is described here by 



I(r) = - 
n 



i-ri 



(5) 



where l{f) is normalized to unit total flux. The relation between 
T and the more commonly used parameter a, defined by 



(6) 



7(r)oc l-a(l- Vl -r 2 ), 

is given by a = 

The parameterization involves basic microlensing param- 
eters: to (time of closest approach), u a (impact parameter), ?e 
(crossing time of the Einstein ring radius), the source size p* 
(in units of 6b) and two annual parallax parameters tte and iff 
(see Sect. [3 common to all data sets, and specific parameters 
for each site: baseline and blending magnitudes, as well as the 
linear limb-darkening coefficients F. The values found for the 
latter are more fully discussed in Sec. [6] Errors on the best- 
fit parameters were obtained by Monte-Carlo simulations: we 
generate 500 sets of noisy light curves for each of the five ob- 
serving sites using the obtained best-fit parameters, based on 



the rescaled error bars of the data; then, from the distribution 
of the obtained values we determine the (non-Gaussian) 68.3 % 
confidence intervals for each parameter. The best correspond- 
ing set of parameters fitting the data and their errors is given in 
Tab.ffl 

Table 1. Parameters for the best-fitting point-lens limb- 
darkened extended source models, using linear limb-darkening 
law (one LLDC coefficient per data set). The set of data 
contains 612 measurements coming from 4 PLANET sites 
(Danish 1.54m, UTas lm, Boyden 1.5m and SAAO 1.0m) and 
from OGLE. We also report the x 1 value. 



Parameters 


Value & Error 


t (days) 


3166.8194 ±0.0002 


U a 


4.60!™ x 10- 3 


t E (days) 


1 a 9 q +0.04 


P 


4.00+ , , , x 10- 2 


T/Isaao / fl ;lsAAO ■ • • 


35 +o: ° 2 / 45 +a03 

"■-'-'-0.05 ' "■ HJ -0.06 


T/luTas / A/luTas • ■ • • 




F/lBoyden / fl /boyden • 


0.58!°;°°/ 0.68!°;°° 


T/loGLE / a /loGLE ■ ■ • 


n 40 +0.04 , n. c +0.04 
"•^-'-0.06 ' "•-'-'-0.06 


r« banish / a R banish • • 


61 + ° 03 /0 70+ 03 

"• ul -0.06 ' "-'"-0.05 


^b/^sIsaao 


0.25 ± 0.02 


Fb/Fs banish 


0.05 ± 0.01 


^B/^sluTas 


0.55 ± 0.01 


Fb/Fs boyden 


1.15 ± 0.2 


^b/^sIogle 


0.00 + 0.0 


e ii.o.r 


1326/593 



The blending fraction significantly varies from one 
PLANET site to another. This is a consequence of the fact 
that the source has two very close neighbours, one of them at 
a distance of 1.4" as stated in Sec. 12.21 Our measured blend- 
ings are consistent with the different telescopes CCD resolu- 
tions. Moreover, the low value obtained with the Danish 1.54m 
(^ 6%) favors a low blend flux from the lens star. 

Based on this model, we find the transit time to be 
1.05 day and the date (where we write t = HJD - 2,450,000) 
at which the lens starts to transit the source disk, f en tiy - 
3166.2937, as well as the exit time f ex j t 3167.3451. One can 
note the UVES spectra (Sec. 12. 3> were taken after the end of 
the transit, and so only weak differential magnif ication over the 
source disk is expected JHevrovskv et all2000l) . 



4. Nature of the source star 

Measuring the relative proper motion and velocity at the lens 
distance requires a determination of the distance and radius of 
the source star. This can be achieved by determining the appar- 
ent magnitude at baseline, the effective temperature, the grav- 
ity, and the amount of interstellar absorption toward the source. 
We seek to measure these quantities by use of a combination of 
the obtained photometry and spectroscopy. 
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4.1. Source star calibrated colour-magnitude diagram 

Fig. |3 shows the calibrated OGLE (V - 7) versus 7 colour- 
magnitude diagram of the field of the source/lens. The cali- 
brated magnitude and colour of the target are 7 = 16.59 ± 0.05 
and (V - 7) = 2.26 ± 0.08. A box encompasses stars belonging 
to the red giant clump (RC), and the two histograms give a cal- 
ibrated mean position of the clump of 7 RC = 15.80 ±0.10 and 
(V - 7) RC = 2.13 + 0.05. The shift in position of our target with 
respect to the mean red clump position is then AI = 0.79 ±0.11 
and A(V-I) = 0.13 ±0.09. 

For the absolute magnitude of the clump, we adopt th e 
Hipparcos position as given bv IStanek & Garnavichl Jl99fl. 
Mi = -0.23 ± 0.03. The mean Hipparcos clump colour of 
(V - I)o RC = 1-00 ± 0.05 is adopted, assuming any difference 
with the Baade's Window mean colour to be due to the anoma- 
lous extinction la w in the Galactic Bulge, as first explained by 
IPonowskil teOOOl). 

The Galac tic Center distance modulus from 

lEisenhauer et alJ J2005I) is adopted as // = 14.41 ± 0.09, 
but a slightly larger distance is adopted for the red clump in 
this field, due to its negative longitude and the bar geometry 
(see, e.g.: IStanek et all 1 1 9971) . The difference in dista nce is 
based on mean clump position in different regions from lSumil 
and amounts to A/u = 0.21 ± 0.10, giving an adopted 
distance to the red clump of this field of // = 14.62 ± 0.13. 
Therefore, the mean dereddened magnitude of the clump is 
7 Q RC = 14.39 ± 0.14. Comparison with its calibrated apparent 
magnitude and colour gives an estimate of reddening param- 
eters as Ai = 1.41 ± 0.17 and E(V - I) = 1.13 ± 0.07, the 
latter in very good agreement with Sumi's measurement of 
E(V-I) = 1.17 for the nearby OGLE-II field BUL-SC25. From 
these two independently determined values follows the total- 
to-selective absorption ratio Ri = Ai/E(V - 7) = 1.25 ± 0.14, 
intermediate between standard reddening law value (1.5) and 
Sumi's anomalous mean val ue for the Gal a ctic Bu lge (0.964). 

We furthermore use Idirardi et alJ J2002h theoretical 
isochrones to determine the mean metallicity of the red clump, 
and by adopting an age of 10 Gyr (typical age of red clump 
giants in the Bulge) we find Z RC 0.008. 

Finally, from the adopted position of the clump and the shift 
of the target, we obtain a dereddened magnitude and colour of 
the target as 7 = 15.18 ± 0.18 and (V - 7) G = 1.13 ±0.11. 

The comparison of (V - 7 ) with the values from t he grid of 
model atmosphere colours bv lBuser & KuruczH 19921) provides 
an estimate of the photometric temperature of the source star. 
First, we note that typical changes in the theoretical value of 
(V - I) for a change in metallicity by a factor of two, or a 
change in logg by an order of magnitude, are about 5 x 10~ 3 
(for fundamental parameters in the range of values of interest in 
the present study). In comparison, a change in the value of r e ff 
by 100 K gives rise to a 10 times larger change in (V - 7) . The 
value of (V— 7) is therefore primarily a measure of the effective 
temperature of t he star. Interpolat i ng in t he values of (V - 7) D 
from the grid of iBuser & Kuruc3 dl992l) . gives us the best fit 
effective temperature of the source star from the photometry 
alone, as T^ ot = 4500 ± 250 K. The quoted uncertainty reflects 
solely from the uncertainty in the estimate of (V - 7) . 
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Fig. 2. (a) Calibrated (V - 7), 7 colour-magnitude diagram of 
the field around OGLE 2004-BLG-254 as obtained with the 
OGLE telescope at Las Campanas. The target is plotted as a 
red open star. The mean position of the red giant clump is the 
blue dot, as determined by the color histogram (b) and mag- 
nitude histogram (c) of the stars inside the black box of panel 
(a). Mean positions and uncertainties are the blue and red lines 
respectively. 



4.2. Source star fundamental parameters from 
spectroscopy 

The position of the source star (and lens) is in the direction of 
the Sagittarius dwarf galaxy, and it could either be a member 
of the dwarf galaxy or of the Galactic Bulge. The line positions 
of the observed spectrum have a general offset of +134 km/s 
compared to laboratory data, which is also consistent with the 
star being in either the Sagittarius dwarf galaxy or the Galactic 
Bulge. The mean chemical abundance [Fe/H] ranges between 
-0.8 a nd -1.2 in Sagittarius, and around -0.1 in the Galactic 
Bulge jFulbright etal]l200d) . For the analysis of the UVES 
spectrum, we have therefore computed a grid of model atmo- 
spheres with T & ff between 4000 K and 4600 K, log g between 
0.0 and 3.0, and scaled solar abundances with [Fe/H] between 
-2.0 and +0.5, and corresponding synthetic spectra. 
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Fig. 3. A comparison of the observed de-redshifted OGLE 2004-BLG-254 source star spectrum (grey) with a synthetic spectrum 
(thick black line) based on a model atmosphere with T e g = 4100 K, log g — 1 .9 and Z = 2. 1 Z Q . The thin grey line around y = 
indicates the residual. The spectra are normalized to the local continuum. From top to bottom, the panels show the spectral region 
of: (2) a triplet of magnesium lines; (3) a triplet of chromium lines; (4) the NaD doublet; (1) and (5) show two other typical 
regions of the spectrum. 
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The synthetic spectra are calculated by performing 1- 
dimensional radiative transfer through standard Kurucz' model 
atmospheres and disk-integrating the results from 7 yu-angles. 
Spectral line data is taken from the VALD data base 
jKupka et al.lll999T) . The observed spectrum shows approxi- 
mately 10,000 well-defined lines, and almost all of them are 
identifiable from comparison with line positions and strengths 
of the transitions listed in the vald data base. Line profiles are 
computed as Voigt profiles with the necessary broadening pa- 
rameters taken from the data base. 

Among the many atomic lines, we have selected 3 particu- 
larly well suited systems of strong Mg, Cr, and Na lines, whose 
intensity and line shapes are fitted to constrain the possible es- 
timates of the fundamental parameters: effective temperature 
T e ff, gravity logg, and metallicity Z. Other lines are used to 
control this estimate and to get a feeling for individual devia- 
tions from a scaled solar abundance. Fig.[3]shows the observed 
spectrum in those regions, together with a synthetic spectrum 
discussed below. 

Magnesium lines - The triplet of neutral magnesium lines 
around 5175 A is well suited to obtain limits on the tempera- 
ture and gravity. The line system overlaps with the position of 
a relatively strong MgH band, and the ratio between Mg and 
MgH is sensitive to temperature as well as gravity. For strong 
gravity, the atomic magnesium triplet lines become very broad, 
but for low temperatures, the balance shifts in favour of MgH. 
Therefore, the shape and the intensity of the atomic lines can 
be used together with the ratio (or absence) of the intensity of 
MgH relative to the intensity of the atomic Mg lines to provide 
information on both temperature and gravity. The absence of 
MgH in our observed spectrum allows us to conclude that the 
star is not cooler than 4000 K. The breadth of the atomic Mg 
lines allows us to confine the value of log g between 1 .0 and 
2.5. The medium-strong neutral atomic Mg line at 5711 A is 
known to respond oppositely to the triplet lines to changes in 
gravity, i.e. to become stronger for decreasing gravity. The 
synthetic line is too strong for log g = 0.0 and solar metallicity, 
while it becomes too narrow for high log g. 

Chromium lines - As for the Mg triplet region, the MgH 
molecular system also has a relatively strong band in the 
region of a triplet of three strong chromium lines at 5204.51, 
5206.04 and 5208.42 A, which limits r eff to no less than 
4000 K. Models of r eff = 4200 K fit the Cr lines well, while 
models of r e fj = 4400 K result in wings of the lines being too 
weak even for high metallicity models, while = 4000 K 
would require a metallicity considerably below Z Q . The 
chromium system is less sensitive to gravity than the other 
two line systems, and for some T^, even values as low as 
log g ^ 0.0 could be in agreement with the observed spectrum. 
On the other hand the lines are sensitive to the chromium abun- 
dance and Z = 0.3 Z is too low unless is as low as 4000 K. 

Sodium lines - The intensity and form of the NaD lines 
around 5890 - 5896 A and other neutral sodium lines are very 
sensitive to T e ff as well as to gravity and (sodium) abundance. 
Often these lines are not useful for determination of the 
fundamental parameters and abundances, because interstellar 
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Fig. 4. Using x 1 fitting of the spectrum to determine the most 
likely values of Z, T e g and log g reveals a strong correlation be- 
tween these parameters. Here we show the location of the min- 
ima in the;^ 2 surface for five different chosen values of [Fe/H]. 
Around each minimum we indicate the contours correspond- 
ing to confidence intervals of 68.3%, 95.4% and 99.7% (1, 2 
and 3 cr). Treating Z as a free parameter, we found the best^ 2 
solution around T e s = 4100 K, \ogg =1.9 and Z = 2.1 Z Q . 



absorption saturates or changes their intensity. In this case, 
however, the main component of the interstellar absorption is 
redshifted by a velocity of 122 km/s relative to the star, and 
the intrinsic stellar NaD lines are very strong and appear to be 
only moderately affected by interstellar absorption. The fact 
that the fundamental parameters derived from the NaD lines 
are in good agreement with the parameters derived from the 
other stellar lines, also indicates a small interstellar absorption 
at the radial velocity of the star. Model spectra from our grid 
with high metallicity (Z = 3 Zq), high gravity (log g = 3) or 
low T e ff, all give far too broad lines compared to the observed 
spectrum, and can therefore be excluded. Models of low 
gravity (~ 0), high T eS (~ 4400 K) or low Z (~ 0.3 Z Q ) give 
on the other hand too narrow lines, and would require a strong 
interstellar component at the same radial velocity as the star. 

The intensity and form of the line systems discussed above 
are all different functions of T e g , log g and Z, and in princi- 
ple three systems (like the Mg, Cr, and Na systems) are suf- 
ficient to determine the three fundamental stellar parameters 
uniquely from the observed spectrum. In practice, of course, 
the observed spectrum is noisy, and the theoretical spectrum 
suffers from inaccuracy in the model structure, incompleteness 
of the line list, etc. Therefore, rather than a unique fit, there 
is a certain range of models which give acceptable fits to the 
spectrum. If T e ff = 4300 K is adopted, the Mg lines are well 
reproduced by log g — 2 and Z = Z , while the Cr and Na lines 
and most of the other lines are slightly too weak in the wings. 
The fit to the Cr and Na lines can be improved by decreasing 
the effective temperature to 4200 K, although this slightly de- 
creases the goodness of the fit to the Mg lines. A similar effect 
can be obtained by keeping the value r e jf = 4300 K, but in- 
creasing log g a bit, for example to log g = 2.25, and increasing 
the metallicity. At r eff = 4400 K, the fit to the Mg lines would 
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still be correct, but the fit to the Cr and Na lines would be worse 
than for = 4300 K, and a compensation by increasing the 
gravity or the metallicity would require larger values than for 
T e s = 4300 K, and therefore result in a larger decrease in the 
goodness of the fit to the Mg system than for the r e ff = 4300 K 
model. We therefore conclude that there is no consistent solu- 
tion to the fit for T e s = 4400 K, and this value is therefore 
too large, independently of the adopted values of log g and Z. 
On the other hand, changing the temperature in the other direc- 
tion, for example to T e s = 4000 K, a good fit to the spectrum 
would require that the value of the gravity and the metallicity 
be decreased. 

We conclude that good fits to the Mg, Cr, and Na lines (and 
other features in the spectrum) can be obtained for models with 
r eff ranging from 4000 K to 4300 K, log g from 1 .5 to 2.5 and Z 
from 1 .5 Z Q to 2.5 Z , which classify the star as a normal K3 III 
red giant. We determine the parameter values that agree best 
with the observed spectrum by performing a multi-dimensional 
^-minimization using our grid of synthetic spectra covering 
the range of possible valus of r e ff , log g and Z. The optimal so- 
lution has^ 2 /1504 = 2.997, T eff = 4100 K, logg = 1.9, and 
Z = 2.1 Z Q , with a covariance of T e g/ \ogg = 0.787 and cor- 
relation coefficients of T e ff/\ogg = 0.870, Z/logg = 0.679, 
and Z/r e fr = 0.811. Hence, we further conclude that there is 
a strong correlation between T e « , log g and Z in such a way 
that increasing T e ff must be followed by a increased value of 
log g and Z and vice versa. A synthetic spectrum with the opti- 
mal parameters is shown, together with the observed spectrum, 
in Fig. [3] In addition, we show as an illustration of the strong 
correlation between parameters in Fig. @]a contour plot of the 
^-surface for different values of Z in the plane of (T e e, \ogg). 

4.3. Source radius and distance by combining 
spectroscopy and photometry 

In general terms, photometry has a major advantage over spec- 
troscopy in a higher obtainable accuracy of the broadband 
colours than what is obtainable from integrating the spectrum 
(whose slope is often ill defined). With good transformation 
between colours and stellar effective temperatures, an accurate 
value of T e ff is therefore obtainable. However, a major diffi- 
culty with photometry of distant stars in the Galactic plane 
is a large uncertainty in the estimated reddening. We have 
partly got around this problem by applying the red giant clump 
method, which reduces the reddening problem to a question 
of the reddening of the target relative to the red giant clump, 
plus a problem of defining the clump position in an observed 
colour-magnitude diagram of the field around the target. 

The spectroscopic method, on the other hand, has the ad- 
vantage of being able to simultaneously estimate the stellar ef- 
fective temperature and gravity (and metallicity). The main dis- 
advantage, for the present purpose, is an ambiguity (i.e., corre- 
lation) between the obtained values of T e ff and log g. We have 
illustrated this correlation in Fig. |4] for the present spectrum. 
By demanding that our estimated stellar effective temperature 
satisfy both photometry and spectroscopy, we take advantage 
of more information than from one of the methods alone. The 



value of log g can then be restricted further than was possible 
from spectroscopy alone, by restricting the estimate of log g to 
the part of Fig.|4]that overlaps with the photometric estimate of 

Teff- 

By combining the previous photometry and spectroscopy 
analyses, we get a very small overlap between the two ranges 
of effective temperature values, namely r e ff in the interval 
4250 - 4300 K, although it is probably realistic to enlarge this 
estimated interval a bit. 

4.3.1 . Source distance 

The combination of the fundam ental source star pa rameters 
with the theoretical isochrones of iGirardi et al.l d2002l) and our 
calibrated colour-magnitude diagram can provide a measure of 
the source star distance Ds. For a g iven range of T e ff, log g and 
Z, we get from lGirardi et al.l (120021) isochrones the correspond- 
ing values of the source absolute magnitude Mj\ the distance 
then comes from Mj — I — fi, where fi — 5 log Ds/10 pc is the 
distance modulus. We also require the source to be part of the 
Bulge, for the probability of observing a microlensing event for 
a source outside of the Bulge is very low. 

By assuming a temperature in the range T e ff = 4250 - 
4300 K, we find no solution, whatever the age, gravity or 
metallicity, as illustrated in panel (a) of Fig. [5] This discrep- 
ancy between photometric and spectrocopic temperature can- 
not originate from a higher magnification of the limb of the 
source, since the transit of the lens was already finished when 
the UVES spetra were obtained and the corresponding effect is 
weak (Sec. iFulbright et all feOOdl also reported discrep- 
ancies for Bulge K-giants, with systematic differences in pho- 
tometric and spectrocopic temperatures. 

We then examine the case of higher temperatures (T e ff > 
4300 K). We first remark that the larger the distance, the 
smaller the temperature has to be. the optimal solution is then 
found for the maximum possible distance compatible with the 
source beeing a Bulge star, which we assume to be Ds = 
10.5 kpc. Moreover, this compromise is better w hen the age of 
the sou rce is maximal, which means 16 Gyr in lGirardi et alJ 
d2002l) isochrones. A possible blend of the source star by the 
lens would also imply even more discrepant results, and fol- 
lowing Sec. 13.21 we choose to neglect it. We study now the re- 
sults involving different metallicities. By setting for the source 
the red clump metallicity, Z 0.008 (Sec. 14. 1> . we obtain a 
solution with T e ff 4575 K and \ogg 2.36, and a source 
at the distance Ds - 10.5 kpc. Although T e ff is still com- 
patible with our photometric analysis of Sec. 14.11 it is above 
the higher bound of the allowed spectroscopic range of tem- 
peratures. The corresponding solution is shown in panel (b) of 
Fig. If we adopt a solar metallicity, then our requirement to 
have Z?s < 10.5 kpc leads to a best solution with T e s - 4620 K, 
logg =i 2.31 withDs = 10.5 kpc, which implies similar conclu- 
sions than when assuming the red clump metallicity. This case 
is shown in panel (c) of Fig. [5] We do not find any satisfying 
solution for higher metallicity. 

In conclusion, we adopt the maximum distance compatible 
with the source star being part of the Bulge, Ds - 10.5 kpc, 
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as a best compromise between photometric and spectroscopic 
analysis. 

4.3.2. Physical radius of the source 

The value of the source distance Ds derived above combined 
with the angular size of the source 0» is now used to determine 
the physical radius of the source, fi, = 6, x Ds. 

Th e surface brightness relation from iKervella et alJ 
II20041) 1 : 



log(0,///as) = 3.212 - 0.2 I a + 0.421 (V - /) 



(7) 



links the apparent angular radius of the source to the measured 
dereddened magnitude and colour of the source star. Adopting 
the values of Sec. 14.11 we find an apparent angular radius of 
0» = 4.5 ± 0.7 //as. 

If we include a possible difference in reddening between 
the source and the red giant clump, we can write the real dered- 
dened magnitude and colour of the source as: 



h = 15.18 ± 0.18 -RiAE(V-I), 
(V-I) a = 1.13 ±0.11 - AE(V- 1), 



(8) 
(9) 



where A is in the sense source minus red clump. From this, we 
get an apparent angular radius: 



log(0,//ias) = 0.652 -kAE(V-I), 



(10) 



where the coefficient of the differential extinction Ic = 0.421 - 
0.2 Ri reads 0. 17 + 0.05 for our adopted value of Rj. However, a 
differential reddening is unlikely, since it would imply a smaller 
reddening for a star behind the red clump, which is somehow 
unsatisfactory. 

Thus, assuming no differential extinction, we derive a value 
of R t =i 10.2 R . We also checked that by using the simple 
relation L oc R^ x T 4 , where T and L are deduced from the 
isochrones for a source star fulfilling the conditions shown in 
panel (b) or (c) of Fig. [5] we find a similar result (±2%) for the 
source radius. 

5. The lens location 

With the angular Einstein radius being related to the angular 
source radius 0» as 0e = 0*/p*, we find 0e - 1 12//as. This en- 
ables us to calculate the the relative lens-source angular proper 
motion: ft = 0e/?e - 3.1 mas/yr. The relative velocity v be- 
tween lens and source at the lens distance then follows from 
v = D L /i. 

In principle, a measurement of the source size both in 
Einstein radius and physical units, as well as the measurement 
of parallax parameters completely determine the lens location 
(given the source distnce Ds). 

1 Although this calibration concerns Cepheids, it has been repeat- 
edly demonstrated that surface-brightness relations f or stable giants 
and Cepheids agree within 1 % jNordsren et aj2002h . A possible al- 
ternative would be to make use of surface-brightness relations directly 
calibrated for gian t stars, but they only exist in (V - K). Such a re- 
cent calibration bv lGroenewegenl 120041) leads to an angular radius of 
4.4 yuas, in good agreement with our adopted value. 
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Fig. 5. Same colour-magnitude diagrams as in Fig.|2] The de- 
reddened red clump is the red triangle, the source is the blue 
dot. Isochrones with metallicity Z given in the right corner of 
each figure are plotted in light orange. Photometric temperature 
J e ff and source distance are also indicated. An agreement be- 
tween photometric and spectroscopic observations is obtained 
when the source lies on the isochrone, the vertical position of 
which determines the source distance. In panel (a), it is seen 
that by setting the range of temperatures T e g = 4250 - 4300 K 
from the spectroscopy analysis, no satisfying solution can be 
found. In fact, the two bigger dots on the isochrones match- 
ing this criterion would require a redder source and a distance 
which would put the source outside of the Bulge. Moreover, 
since the absolute magnitude of the target already agues for 
a source on the far end of the Bulge, a differential magnifi- 
cation of the source relative to the clump would require the 
isochrone to be at the left side of the blue dot, which would lead 
to even higher distances. If we allow a higher temperature, so- 
lutions can be found with our requirement that Ds < 10.5 kpc. 
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The source-size parameter p, (in Einstein units) is well- 
constrained by our photometric model (Sec. l3.2t . and from the 
value of deri ved in Se c.EjXH we obtain a constraint on the 
lens mass M as (bominikll998l) : 



— (*) = 



1 Rl 



AGM D S pl 
where x = D^/Ds- 



1 -x 



(11) 




7 8 9 

D L [kpc] 




Fig. 6. Probability densities for the lens mass M, the lens dis- 
tance £>l and the relative transverse velocity v at the lens dis- 
tance, assuming a source distance of Ds = 10.5 kpc and the 
adopted Galaxy model. 



Similarly to Eq. ( fTTl . a measured parallax parameter would 
provide a relation between lens mass and x, but unfortunately 
the event is too short (~ 13 days <k 1 year) to provide a mea- 
surement of the parallax, or even to give reasonable limits: val- 



ues as different as n E ~ 0.01 and 10 are hardly distinguish- 
able from the light curve fit. For a very similar event (duration 
?e - 13 day s, source size p„ ^ 0.06 and very small impact pa- 
rameter mA IYoo et alJ J2004I) only obtained a marginal parallax 
measurement too. 

We then use es timates of the physical parameters, follow- 

[I I I a j a 

Dominik (2006) and assuming his adopted Galaxy model. 

The event time-scale fe = 13.23 days and the source-size pa- 
rameter p„ = 0.04 provide us with probability densities for the 
lens mass M, the lens distance and the relative transverse 
velocity v at the lens distance, which are shown in Fig. [6] 

From these, we find a lens mass M ^ 0.11 M Q , a velocity 
v - 145 kms , and a lens distance Z\ — 9.6 kpc for an as- 
sumed source distance of Ds = 10.5 kpc. The lens is preferred 
to reside in the Bulge, with 58 % probability, even though with 
a source on the far end of the Bulge, disk lenses play a promi- 
nent role with 42 % probabilty for a disk lens scenario. 

Finall y, similarly to the d iscussion of event OGLE 2002- 
BLG-069 JKubas et alJl2005h . we put upper limits on the lens 
mass based on the measured blend ratio F^/Fs, assuming alu- 
minous main-sequence lens star and a source at Ds ^ 10.5 kpc. 
The latter limits, presented in Fig. are compatible with the 
lens being an M dwarf in the mass regime derived in Fig.|6] 
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Fig. 7. Assuming the fraction of blended light being solely due 
to the lens, upper limits on the lens mass can be derived from 
the measured blend ratio Fb/Fs- The implied blend ratio for 
range of main sequence lens stars of spectral types M0-A9 
is plotted as a function of the lens distance (lower horizontal 
axis), or equivalently as a function of the lens mass (upper hor- 
izontal axis) given Eq. (^3. The derived values are compatible 
with the M dwarf lens of mass ~ 0.05 - 0.3 M Q . 



6. Limb-darkening measurements 

The photometric data of OGLE 2004-BLG-254 were dense 
enough to measure limb darkening of the source. Because 
of the duration of the event, the data coverage from a sin- 
gle site is not sufficient for such a measurement. Here, we 
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Fig. 8. Linear limb-darkening coefficients (LLDC) for the I- 
band data sets (upper panel, from bottom to top: SAAO 1.0m, 
OGLE, UTas lm and Boy den 1.5m) and the R-band data set 
(lower panel, Danish 1.54m). The colors follow the same con- 
vention as for Fig. ^ The best value for each r-coefficient 
is indicated by an horizontal line, where the filled area en- 
close the one-cr error bar for each parameter. For the I-band, 
UTas lm and OGLE give consistent values, while SAAO 1.0m 
and Boyden 1.5m are well below or above the latter ones 
(see text for a detailed discussion). Model atmosphere LLDCs: 
filled red dots correspond to the Claret (2000) LLDCs while the 
green squares correspond to the "new fit" LLDCs (see text), 
both for a range of temperature, r eff = 4000, 4250, 4500 K 
(in abscissa), logg = 1.5, 2.0, 2.5 and [Fe/H] = 0.0, 0.1, 0.2, 
0.3, 0.5, compatible with OGLE 2004-BLG-254 stellar param- 
eters. 



benefit from our round-the-clock follow-up that permits mon- 
itoring the event over the full course of the caustic pas- 
sage. We find that a square root limb-darkening law I(r) = 
i[l-r(l- | Vl -r 2 )- A(l - f (1 -r 2 ) 1/4 )] does not im- 



prove the fit, and the strong correlation between T and A leads 
to an unsatisfactory and ambiguous result. 

The limb-darkening coefficients (hereafter LLDC) derived 
from our best fitting model for the I- and R-band are presented 
in Tab. [Hand Fig. [8] For each data set, the best fitting value is 
plotted as an horizontal line, within the a filled area that encom- 
passes the allowed values at the one-cr error level. The color 
conventions are identical to those of Fig. ^ We first note the 
very good agreement between the I-band OGLE and UTas lm 
LLDC measurements, while SAAO 1.0m and Boyden 1.5m 
are not consistent with the latter values. We justify this by 
the fact that the peak region best constrains the LLDC, es- 
pecially the shape of the top of the light curve. As a matter 
of fact, UTas lm shows a good coverage of this region, as 
well as OGLE and Danish 1.54m, as seen in Fig.[0 However, 
even though SAAO 1.0m should also give good constraints, a 
high scatter in the peak region ruins the reliability of the mea- 
surement, and Boyden 1.5m exhibits large scater and has no 
baseline data points. Thus we posit the LLDC coming from 
UTas lm and OGLE (for the I-band) and Danish 1.54m (R- 
band) to be precise enough to be compared with atmosphere 
model LLDC. As a result, the LLDC values we find for the 
two considered bands are in disagreement with the atmosphere 
model. The discrepancy is reduced when higher temperatures 
are considered, which tends to reinforce the conclusion from 
stellar parameters analysis in Sec.|4] 

Fig. [8] then illustrates the comparison between LLDC 
ATLAS atmosphere models predictions with OGLE 2004- 
BLG-254 measurements. The filled red dots in Fig. [8] corre- 
spond to the LLDCs predicted by Claret (2000), for a range 
of paramater values, i.e. r eff = 4000, 4250, 4500 K (in ab- 
scissa), logg= 1.5, 2.0, 2.5 and [Fe/H] =0.0,0.1,0.2,0.3,0.5, 
compatible with the stellar parameters of OGLE 2004-BLG- 
254. However, we point out that Claret's values are obtained 
by least-squares fits of intensity points more or less regularly 
spaced in incidence angle fi = Vl - r 2 . In terms of the radial 
position on the stellar disk r, such a fit gives very high weight to 
points close to the limb and very low weight to points close to 
the center. In order to avoid this bias, we interpolate the ATLAS 
points by cubic splines and perf orm the least-squar es fits on the 
obtained I(r) intensity profiles jHevrovskvll200l . These new 
LLDC values are plotted as green squares in Fig. [8] As shown 
in the last column of Tab.|2]for the given sources, the resulting 
LLDCs are systematically lower than Claret's, in this parame- 
ter range by about 0.05. As seen from Fig- (HI sucn a difference 
in the LLDC may correspond for example to a temperature dif- 
ference of several hundred K. 

Previous microlensing events have yielded nine limb- 
darkening m easurements by sever a l follow-up teams: six for 
GK giants J Albrow et alJ Il999bl l200ot iFields et alJ 120031: 
ICassan et al.l2004Hjiang et alj2004HYoo et alj2004 . one for a 
subgiant jAlbrow et alj20 0l|). and tw o for main sequence stars 
dAfonso et al.l2000HAbe et alJ2003h . The /-band GK giant re- 
sults are the most relevant for comparison, and are summarized 
in Tab. giving the r eported LLD C alongside values derived 
from ATLAS models dClareJ2000l) for comparison. 
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Event 



Source characteristics 
Type T eff logg 



[Fe/H] 



Measured LLDC 
ai 



ATLAS LLDC 

a. 



ATLAS 
a/ (new fit) 



MACHO 1997-BLG-28 
MACHO 1997-BLG-41 
EROS 2000-BLG-5 
OGLE 2002-BLG-069 
OGLE 2003-BLG-262 
OGLE 2003-BLG-238 
OGLE 2004-BLG-254 



K2 III 
G5-8 III 
K3 III 
G5III 
Kl-2 III 
K2 III 
K3 III 



4250 K 
5000 K 
4200 K 
5000 K 
4500 K 
4400 K 
4000 - 4300 K 



2.0 
3.2 
2.3 
2.5 
2.0 
2.0 
1.5-2.5 



0.0 
-0.2 
+0.3 
-0.3 

0.0 

0.0 
+0.2-+0.4 



0.83 (±15%) 
0.46 
0.54 
0.60 
0.70 + 0.13 
0.57 ± 0.06 
53 +0 04 

-0.06 



0.65 
0.58 
0.67 
0.57 
0.63 
0.64 
0.63 (closer) 



0.60 
0.53 
0.62 
0.52 
0.58 
0.59 
0.58 (closer) 



Table 2. Limb-darkening coefficients for the /-band of GK Bulge giants in OGLE 2004-BLG-254 and published events. The 
column ATLAS LLDC is from Claret (2000), whereas the last column is a new fit to ATLAS model atmosphere intensities. The 
best measured LLDC for OGLE 2004-BLG-254 correspond to 7/ eff = 4500 K, logg = 1.5 - 2.5 and [Fe/H]= 0. - 0.1, based 
on the LLDC derived from UTas and OGLE data (see text). The ATLAS values given here are the closest to the measured ones 
(cf. Fig.|8}. 



We exclude three cases from Tab. [2] for our comparison: 
OGLE 2002-BLG-069 because it is of earlier spectral type; 
MACHO 1997-BLG-41, also of earlier type, for it suffers 
from correlations between the LLDC and other model pa- 
rameters, and OGLE 2002-BLG-262 which has large uncer- 
tainties due to sparse temporal coverage. The four remain- 
ing objects are all K giant stars. The comparison between 
LLDC measurements and ATLAS atmosphere models are pre- 
sented in Tab. |2] MACHO 1997-BLG-28, involved a cusp 
passage and seems to deviate from theory. While we do not 
mistrust the underlying light curve model, the accuracy of its 
derived limb-darkening coefficient is not sufficient to chal- 
lenge the atmosphere modelling: an uncertainty of 15-20 % 
seems reasonable and removes part of the apparent discrep- 
ancy. Assuming Claret (2000) limb-darkening coefficients for 
the stars with their known temperature, the two remaining liter- 
ature events, EROS 2000-BLG-5 and OGLE 2003-BLG-238 
- together with OGLE 2004-BLG-254 -, appear to disagree 
with theory predictions, considering the high quality and suit- 
ability of the data for a direct surface-brightness profile con- 
straint. However, assuming our "new fit" limb-darkening co- 
efficients, OGLE 2003-BLG-238 now agrees with the pre- 
diction, EROS 2000-BLG-5 is much closer to the predicted 
model, and the discrepancy for OGLE 2004-BLG-254 is re- 
duced. 

There is still an incompatibility between the measured and 
theoretical limb darkening of Bulge K-giant sta rs, as previ- 
ously for EROS 2000-BLG-5 JFields et alJ2003h . However, if 
a limb-darkening coefficient following our fitting prescription 
is used, the discrepancy is significantly smaller than with the 
coefficients published by Claret (2000). 

Finally, a comparison between the shape of synthetic at- 
mosphere models and linear and square-root limb-darkening 
law curves suggests that the classical laws are too restrictive 
to fit well the microle nsing observations, as already suggested 
evrovs kl il2003l) . For example, all the normalized curves 
derived from the classical laws are constrained to pass through 
the point at r — Y5/3, whereas curves derived from synthetic 
spectra of giant stars tend to intersect ~ 5 % closer to the cen- 



ter. If the latter is closer to reality, the LDC inferred from a 
given microlensing event might be biased by the attempt to 
compensate for the too steep outer behaviour of the linear limb- 
darkening law, and this would only be exacerbated by any spar- 
sity in the photometri c coverage at that phase. We note here 
that IHevrovskvl d2003l) using simulated single-lens microlens- 
ing light curves showed that the accuracy of recovering linear 
limb-darkening coefficien ts is limited by the inadequacy of the 
linear limb-darkening law. IHevrovskvl J2003I) suggests address- 
ing this problem using a principal component analysis (PCA) 
approach, where orthogonal basis functions extracted from a 
grid of atmosphere models are used to describe the broad-band 
limb darkening of stars. This will be investigated in a forthcom- 
ing paper together with a set of similar single lens events with 
finite source effects. 

7. Summary and conclusion 

We have performed dense photometric monitoring of the mi- 
crolensing event OGLE 2004-BLG-254, a relatively short du- 
ration and small impact parameter microlensing event gener- 
ated by a point-like lens transiting a giant star. The peak magni- 
fication was about A ~ 60, effectively multiplying the diame- 
ters of our network telescopes by a factor ~ 8. High-resolution 
spectra were taken using the UVES spectrograph at La Silla 
while the source was magnified by a factor A ~ 20, just af- 
ter the end of the transit of the source over the caustic. This 
yielded a precise measurement of the characteristics of the 
star, a K3 III giant, although a discrepancy is found between 
photometric and spectroscopic temperature. Using a calibrated 
colour-magnitude diagram analysis and isochrones, we find a 
source angular radius 0* 4.5 + 0.7 //as, and a physical radius 
R* lO.2/?0 with the source distance Ds = 10.5 kpc. A Galaxy 
model together with the event time-scale ?e = 13.23 days 
and the source-size parameter p» = 0.04 yielded a lens mass 
M =i 0.11 M Q , a lens distance = 9.6 kpc, and a velocity 
v - 145 km s at the lens distance. From our photometric data, 
we have derived measurements of the source limb-darkening 
coefficients for the / and R broadband filters, and provided ar- 
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guments for a discussion about a lack of relevant physics in 
K-giants atmosphere models, or an inadequacy of the linear 
limb-darkening law to model. 
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